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Abstract 

A  spatial  and  temporal  numerical  simulation  has  been  carried  out  of  a  pulsed 
(100%  modulated),  rf  inductively  coupled  plasma  discharge  in  argon, 
connected  to  an  additional  (ballast)  diffusion  chamber  of  much  larger 
volume.  It  is  demonstrated  that  during  the  active  phase,  the  presence  of  the 
large  ballast  volume  has  a  small  impact  on  the  parameters  of  the  plasma  in 
the  smaller  discharge  chamber.  In  this  case  the  plasma  parameters  in  the 
discharge  chamber  can  be  estimated  separately  from  the  diffusion  chamber 
by  a  standard  method  using  the  characteristic  ambipolar  diffusion  time  (for 
example,  using  a  global  model).  However,  during  the  afterglow  phase,  the 
situation  is  changed  significantly.  In  the  afterglow,  the  densities  of  charged 
particles  in  the  discharge  chamber  become  lower  than  in  the  large  ballast 
chamber  due  to  more  rapid  diffusion  loss.  As  a  result,  the  reverse  of  the 
active  phase  situation  occurs,  namely,  the  plasma  does  not  flow  from  the 
small  to  the  large  chamber,  but  in  the  opposite  direction,  from  diffusive  to 
discharge  volume,  and  both  the  plasma  density  gradient  and  the 
self-consistent  ambipolar  electric  field  in  the  small  chamber  change 
directions.  This  phenomenon  leads  to  new  effects  in  the  discharge  volume, 
in  particular  a  decreasing  rate  of  decay  of  densities  of  charged  particles  and 
electron  temperature.  Thus,  in  the  afterglow  the  presence  of  a  large 
additional  ballast  volume  has  a  significant  impact  on  the  plasma  transport. 

In  this  case,  a  simple  treatment  of  the  plasma  in  the  discharge  chamber  in  the 
framework  of  a  spatially  averaged  model  (for  example,  the  global  model)  is 
inadequate. 


AQ1 


(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

Utilization  of  a  pulsed  or  modulated  discharge  is  a  possible 
approach  for  optimizing  and  improving  the  various  plasmo- 
chemical  reactors  and  devices  for  plasma  processing,  etching, 
production  of  nanotubes,  etc.  By  varying  the  duration  of 
the  active  phase,  repetition  frequency,  deposition  power  and 
other  external  parameters  it  is  possible  to  gain  increased 
control  of  various  plasma  characteristics  and  therefore  achieve 
improvement  in  the  operation  of  these  devices  [1,2]. 


For  describing  and  optimizing  these  systems  both  large- 
scale  computational  codes  (see,  for  example,  [3-12])  and 
simple  analytical  models  are  widely  used.  These  allow 
quick  estimation  of  plasma  parameters,  prediction  of  the 
discharge  behavior  and  developing  scaling  laws  over  a  range 
of  conditions.  In  the  last  case,  the  most  common  approach 
is  based  on  the  utilization  of  spatially  averaged  models  (see, 
for  example,  [33])  when  the  diffusion  of  particles  of  type  j  is 
estimated  through  their  effective  diffusion  time  Xj  —  A2/ Dj, 
where  A  is  the  diffusion  length  of  the  discharge  volume, 
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Figure  1.  Schematic  diagram  of  the  experimental  setup. 

which  for  a  cylinder  is  R/2A,  where  R  is  the  radius  of  the 
cylinder  and  Dj  is  the  diffusion  coefficient  of  the  particle.  For 
more  complicated  geometries,  the  reactor’s  effective  diffusion 
length,  in  the  simplest  approximation,  can  be  estimated  as  a 
ratio  of  the  volume  to  the  surface  area,  A  ~  V/S.  For  steady- 
state  as  well  as  for  pulse  discharges  many  calculations  using 
this  approach  have  been  published,  and  following  a  series  of 
papers,  [14-17]  is  frequently  referred  to  as  the  ‘Global  model’ . 

Application  of  the  Global  model  for  pulsed  discharges 
may  be  somewhat  more  difficult  than  steady-state  situations 
and  it  may  be  important  to  take  into  account  design  features  of 
the  reactor  which,  at  first  glance,  may  not  seem  necessary. 
An  example  of  this  situation  is  described  and  presented  in 
this  paper.  Many  plasma  processing  systems  use  a  relatively 
small  discharge  volume  connected  to  a  much  larger  diffusion, 
or  ballast,  chamber.  In  non- stationary  discharges  the  presence 
of  the  two  significantly  different  discharge  and  ballast  volumes 
may,  under  certain  circumstances,  lead  to  unexpected  results. 
Specifically,  the  ballast  volume,  which  accumulates  energy  and 
particles  from  the  discharge  chamber  during  the  active  phase, 
can  return  these  during  the  afterglow.  The  application  of  simple 
average  models,  which  do  not  take  this  into  account,  can  lead 
to  significant  errors. 

Thus,  the  aim  of  this  paper  is  to  study  the  plasma 
characteristic  in  a  pulsed  rf  inductively  coupled  plasma  (ICP) 
discharge,  initiated  in  a  small  chamber,  which  is  connected  to 
a  diffusion  chamber  of  a  much  greater  volume. 

2.  Experimental  setup  for  simulations 

The  simulations  in  this  work  were  performed  for  the  system 
used  in  our  previous  work  [18-20],  where  details  of  the 
experiment  can  be  found.  A  simplified  schematic  diagram 
of  the  setup  is  shown  in  figure  1.  The  system  volume  is 
divided  into  two  chambers.  The  rf  ICP  discharge  is  created  in 
a  small  chamber  (discharge  chamber  in  figure  1)  with  a  length 
of  8.3  cm  and  a  radius  of  5.5  cm.  The  plasma  can  diffuse  into 
the  large  volume  (diffusion  chamber),  with  a  length  of  122  cm 
and  a  radius  of  30  cm.  Working  gas  is  argon  at  pressures 
between  5  and  30mTorr.  Power  deposited  in  the  discharge 
varies,  typically,  over  the  range  from  5  to  300  W  with  an  active 
phase  of  the  discharge  from  100  to  300  fis  and  an  afterglow 
phase  from  1  to  3  ms.  Current  and  voltage  probes  are  used  to 
measure  rf  power  into  the  coil. 

It  should  be  noted  that  both  probe  (see,  for  example  [21]) 
as  well  as  optical  (see,  for  example,  [22])  diagnostics  are 
well  developed  for  rf  ICP  discharges.  In  particular,  emission 


spectroscopy  allows  spatial  and  temporal  measurement  of 
excited  states  of  atoms  [23,24].  However,  for  ICP  discharges 
of  the  specific  geometry,  including  additional  ballast  volume, 
shown  in  figure  1,  there  are  no  detailed  results  from  either 
experiments  or  simulations.  This  paper  presents  results  of 
simulations  only.  Experimental  confirmation  of  the  effects 
outlined  here  will  be  provided  in  a  future  publication. 

3.  Simulation  model 

A  full-scale  2D  simulation  of  the  above  plasma  reactor  in 
argon  has  been  performed  using  the  methodology  previously 
developed  [25] .  The  main  input  parameters  for  the  simulations 
are  as  follows:  geometry  of  the  discharging  volume,  pressure 
and  composition  of  gas  and  rf  power  deposited  in  the  plasma 
(or  current  in  the  coils).  Early  simulations  showed  that  the 
frequently  used  three  level  scheme,  including  only  the  ground 
state,  one  excited  state,  namely,  the  metastable  level  (see, 
for  example,  [7]),  and  the  ion  ground  state  resulted  in  large 
discrepancies  when  compared  with  measurements  of  afterglow 
parameters.  By  including  a  resonant  level,  energetically  close 
to  the  metastable  level,  acceptable  accuracy  in  the  simulation 
was  achieved.  This  resonant  level  mixes  effectively  with  the 
metastable  level,  which  is  especially  important  in  describing 
the  afterglow.  Therefore,  for  the  simulations  reported  here,  this 
four  level  model  was  used.  The  balance  between  metastable 
and  resonance  atoms  strongly  depends  on  both  direct  excitation 
from  the  ground  state,  but  also  on  cascade  contributions  from 
higher  excited  states.  Cascade  contributions  were  taken  into 
account  using  the  approach  of  Pitchford  et  al  [27],  i.e.  the 
total  flux  from  higher  levels  is  partitioned  into  metastable 
and  resonant  levels  in  a  fixed  ratio  which  can  be  set  in  the 
simulation  (in  the  simplest  case  1:1).  Rate  constants  for 
processes  involving  electrons  and  their  transport  coefficients 
and  in  equations  for  densities  ( ne )  and  average  electron  energy 
(temperature  Te),  as  well  as  densities  of  metastable  ( nm )  and 
resonance  ( nr )  levels,  were  found  using  calculated  electron 
energy  distributions  (EED).  The  EED  has  been  found  from 
the  Boltzman  transport  equation,  taking  into  account  electron 
heating  by  the  electric  field  and  energy  exchange  through 
elastic,  electron-electron  and  inelastic  collisions.  The  set  of 
physicochemical  reactions  used  in  the  simulations  are  shown 
in  table  1 . 

The  set  of  the  ion  balance  equations,  written  in  the  drift- 
diffusion  approximation,  and  the  EED  were  solved  numerically 
using  the  CFDRC  software  [4].  This  code  is  capable  of 
handling  a  variety  of  plasma  systems,  one  of  which  is  the 
pulsed  ICP  plasma.  A  detailed  self-consistent  model  of  the 
discharge  plasma,  numerical  iteration  scheme  and  technique 
for  solving  the  set  of  equations  are  described  in  detail  in  [28], 
and  examples  of  how  this  scheme  is  used  to  solve  different 
plasma  problems  can  be  found  in  [25, 26, 29, 31, 32]. 

4.  Simulation  results 

Typical  results  of  the  2D  simulations  of  the  spatial  distribution 
of  electron  density  at  the  end  of  the  active  discharge  phase 
and  500  jus  in  the  afterglow  are  presented  in  figures  2  and  3. 
As  might  be  expected,  the  peak  electron  density  during  the 
active  phase  occurs  in  the  discharge  (small)  chamber  not  too 
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ICP  discharge  with  a  ballast  volume 


Table  1.  Set  of  plasma  reactions  used  in  the  model. 


#  Fraction 


e  (eV)  Constant 


Comment 


1 

e  +  Ar  -■ 

►  e  +  Ar 

— 

2 

e  +  Ar  ->  e  +  Ar* 

11.5 

3 

e  +  Ar  -■ 

►  e  +  Ar* 

11.67 

4 

e  +  Ar  2e  +  Ar+ 

15.9 

5 

e  +  Arm  ■ 

->  2e  +  Ar+ 

4.35 

6 

e  +  Ar;. 

2e  +  Ar* 

0.07 

7  2Ar^  ->  efast  +  Ar  +  Ar+  — 

8  2Ar*  ->  efast  +  Ar  +  Ar+  — 

9  Ar^  +  Ar*  -*  efast  +  Ar  +  Ar+  — 

10  Ar*  — »  Ar  +  hv  — 

11  e  +  Ar^  —>  efast  — 


Cross  section  [33] 

Effective  cross  section  compiled  from  [33, 34] 
Effective  cross  section  compiled  from  [33, 34] 
Cross  section  [33] 

Cross  section  [36] 

kv  =  1  x  1013  m3  s-1  [35] 

kv  =  1  x  1013  m3  s-1  [35] 
kp  =  1  x  1013  m3  s-1  [35] 
kp  =  1  x  1013  m3  s-1  [35] 

Ar  =  3  x  106  s-1 

Kp  =  4  x  10-16  s-1 


Momentum  transfer 
Metastable  state  excitation  (11.55  eV) 
Resonant  state  excitation  (1 1.67  eY) 
Direct  inonization 
Stepwise  ionization 

Mixing  from  metastable  to  resonant  state 
(quenching) 

Penning  ionization 
Penning  ionization 
Penning  ionization 
Resonance  radiation  with  adjustment 
for  self-absorption  (y  =  106.4  nm) 
Superelastic  collision 


1  l  l  l  l  l  l 
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Figure  2.  Spatial  distribution  of  electron  density  at  the  end  of  the 
active  phase. 
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Figure  3.  Spatial  distribution  of  electron  density  at  t  =  500  jlcs  in 
the  afterglow. 

far  from  the  rf  window  and  due  to  ambipolar  diffusion,  extends 
somewhat  into  the  large  ballast  volume  (see  figure  2).  Again, 
as  one  would  expect,  electrons  are  generated  in  the  small 
discharge  chamber  and  act  as  a  source  for  the  ballast  chamber. 
In  the  afterglow,  the  peak  electron  density  moves  outside  the 
discharge  chamber  into  the  ballast  volume  (see  figure  3).  The 
overall  density  has  decreased  in  the  afterglow,  but  the  peak 
value  is  now  located  in  the  ballast  volume.  In  this  case,  the 


Table  2.  Spatially  averaged  plasma  parameters  for  the  discharge 
chamber  ‘open’  and  ‘closed’  geometry. 


End  active  phase 

Ne  ( m- 3) 

Nm 

(m  3) 

Te  (eV) 

Open 

3.2  x  1017 

2.6 

x  1017 

2.4 

Closed 

4.5  x  1017 

2.8 

x  1017 

2.56 

ballast  volume  acts  as  a  source  of  electrons  for  the  discharge 
chamber.  Thus,  during  the  active  phase  the  ballast  volume  has 
only  a  small  effect  on  the  behavior  in  the  discharge  volume,  but 
during  the  afterglow  phase  the  ballast  volume  is  a  significant 
perturbation  on  the  discharge  chamber. 

In  order  to  demonstrate  this  behavior  more  clearly,  special 
simulations  were  performed  for  the  device  as  shown  in  figure  1 
(henceforth  referred  to  as  ‘open’  geometry)  and  for  the  case 
when  the  opening  of  the  discharge  chamber  is  closed  by  an 
additional  wall  (henceforth  referred  to  as  ‘closed’  geometry), 
for  the  same  excitation  conditions.  Table  2  shows  the  Spatially 
averaged  parameters  of  the  plasma  in  the  small  chamber  at 
the  end  of  the  active  phase  for  the  two  cases.  As  seen  in  the 
table,  the  values  calculated  for  the  open  and  closed  discharge 
chamber,  during  the  active  phase,  are  reasonably  close  in 
value.  Therefore,  during  the  active  phase  the  presence  of  the 
ballast  chamber  has  only  a  small  effect  on  the  parameters  of 
the  plasma  in  the  discharge  chamber.  This  can  be  explained 
on  the  basis  of  the  simple  balance  of  quasistationary  plasma 
during  the  active  phase,  that  is,  the  balance  of  creation  and 
diffusion  decay  of  charged  particles  to  the  walls  (Vi ra  =  1). 
Since  all  of  the  applied  power  goes  into  the  small  chamber, 
even  a  much  slower  diffusion  loss  of  plasma  in  the  ballast 
chamber  cannot  compensate  for  the  low  rate  of  its  production 
(the  only  charged-particle  source  is  the  ‘outflow’  from  the 
discharge  chamber).  Under  conditions  of  equal  applied  power, 
the  differences  between  diffusion  lengths  (we  use  ‘O’  and  ‘C’ 
to  designate  open  and  closed  geometry)  A0  and  Ac  can  be 
approximated  by 


1 


2.4-  t r 


1 


2.42 


— y  =  — —  +  — -  and  — r  =  — —  +  - — (1) 


Ao 


R 2  AH2' 


which  leads  to  Ac  <  Ao  and  their  effective  diffusion  times 
rc  <  To-  The  somewhat  higher  values  for  electron  density 
and  temperature  seen  for  the  closed  geometry  in  table  2  are 
due  to  this  decreased  diffusion  time,  which  for  equal  power 
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results  in  an  increased  Te  and  ionization  rate  (as  the  diffusion 
rate  increases,  the  ionization  rate  must  increase  to  balance  the 
loss). 

Now  consider  the  afterglow  phase,  where,  in  contrast 
to  the  active  phase,  there  is  a  significant  difference  in  the 
plasma  characteristics  with  or  without  the  ballast  volume.  As 
expected,  in  the  simplest  case  (closed),  when  there  is  only  a 
small  chamber  (i.e.  no  ballast  volume)  evolution  of  the  plasma 
density  has  the  well-known  behavior  in  the  form 

d  Ne  Ne  (  t\ 

—  = -  or  Ne(t)  =  Ne( 0)  exp  ( - )  (2) 

At  ra  V  W 

with  the  characteristic  time  for  ambipolar  diffusion  to  the  wall 
(see,  for  example,  [1]) 


with  Da  =  Z>j(l  +  Tq/T\).  This  case  is  well-studied  and  lies 
within  the  framework  of  the  global  model.  A  comparision 
between  results  from  the  simulation  data  and  calculations  using 
equation  (2),  with  temperatures  taken  from  the  simulations, 
agree  very  well.  For  this  simple,  closed  geometry  case, 
calculations  show  the  density  and  temperature  of  the  plasma 
electrons  fall  monotonically,  while  maintaining  a  spatial 
distribution  close  to  the  fundamental  diffusion  mode. 

For  the  case  of  discharge  chamber  with  additional  ballast 
volume,  results  proved  to  be  unexpectedly  different  from 
the  previous,  standard  situation.  In  this  case,  the  spatial 
distribution  of  the  afterglow  plasma  changes  significantly  from 
what  might  be  expected.  Specifically,  the  maximum  in  the 
plasma  density  moves  from  the  discharge  chamber  to  the 
diffusion  chamber,  as  shown  in  figures  2  and  3.  Figures  4 
and  5  show  the  temporal  change  of  average  plasma  parameters 
in  the  discharge  chamber  for  both  the  closed  regime  (dashed 
lines)  and  open  regime  (solid  lines). 

It  can  be  seen  that  the  fall  in  TQ  for  the  closed  geometry 
is  faster  than  for  the  open  case.  As  previously  discussed,  this 
can  be  attributed  to  the  smaller  characteristic  diffusion  length 
for  the  closed  case  and  a  faster  diffusion  cooling  rate.  With 
the  higher  Te  seen  for  the  open  geometry,  one  would  expect 
a  higher  ambipolar  diffusion  rate  and  therefore  a  faster  decay 
of  Ne  in  the  small  chamber.  However,  as  seen  in  figure  5,  the 
decay  of  Ne  in  the  discharge  chamber  is  slower  for  the  open 
geometry  than  the  closed  geometry.  This  indicates  that  the 
discharge  chamber  has  recieved  an  additional  flux  of  charged 
particles  from  the  ballast  volume,  which  compensates  for  their 
decay  at  the  walls. 

Figure  6  shows  the  temporal  (in  the  afterglow)  and  axial 
behavior  of  NQ  in  the  system  for  both  the  open  and  the 
closed  geometries.  The  dashed  lines  correspond  to  the  closed 
geometry,  and  solid  lines  correspond  to  the  open  case.  As 
can  be  seen  from  figure  6,  the  density  of  plasma  in  the  ballast 
volume,  which  is  much  less  than  the  density  in  the  discharge 
chamber  during  the  active  phase,  maintains  this  behavior 
only  at  the  beginning  of  afterglow.  Because  of  a  substantial 
difference  in  the  size  of  volumes,  decay  of  the  plasma  in  the 
ballast  is  much  slower  than  in  the  discharge  volume.  Because 
of  this,  the  electron  density  in  the  discharge  chamber  will  be 
smaller  than  the  density  in  the  ballast  volume  for  only  the  first 
30  /is  into  the  afterglow.  The  plasma  density  gradient  in  the 


Figure  4.  Temporal  behavior  of  averaged  Ne  in  the  discharge 
chamber. 


Figure  5.  Temporal  behavior  of  averaged  Te  in  the  discharge 
chamber. 


Axial  Distance  from  RF  Window  (m) 

Figure  6.  Axial  and  temporal  behavior  of  Nc  for  the  closed  and 
open  geometries.  Times  are  relative  to  rf  cut-off. 

axial  direction  changes  sign  or  diffuses  in  the  opposite  direction 
from  that  seen  in  the  active  phase  (see  figure  7).  The  upshot 
is  that  roughly  speaking,  during  the  active  phase  plasma  flows 
out  of  the  discharge  chamber,  while  in  the  afterglow  it  flows 
back.  Thus,  the  ballast  volume  plays  the  role  of  a  thermostat. 

Lastly,  figure  7  shows  the  temporal  behavior  of  the  total 
ambipolar  flux  crossing  the  imaginary  surface  between  the 
discharge  and  the  ballast  chambers.  Here,  a  positive  flux 
indicates  an  outflow  from  the  small  to  the  large  chamber  and 
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especially  interesting  to  explore  these  phenomena  in  complex 
mixtures,  specifically,  in  electronegative  gases. 
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Figure  7.  Temporal  behavior  of  the  ambipolar  flux  incident  on  the 
imaginary  boundary  between  the  two  chambers. 
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5.  Conclusions 
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pulsed  (100%  modulated)  ICP  in  argon  have  been  performed 
for  the  case  of  a  small  discharge  chamber  connected  to  a  much 
larger  ballast  volume.  During  the  active  phase  of  the  discharge, 
plasma  is  formed  in  the  small  chamber  and  diffuses  to  the 
large  ballast  volume.  Due  to  the  localized  plasma  production 
at  the  rf  window  of  the  small  chamber,  the  plasma  density  in 
the  ballast  volume  is  relatively  low.  As  a  result,  the  presence 
of  the  ballast  volume  is  not  a  significant  perturbation  on  the 
discharge  chamber  plasma  properties  during  the  active  phase. 
This  occurs  despite  a  much  slower  rate  of  diffusion  to  the  walls 
in  the  large  chamber. 

The  situation  is  quite  different  during  the  afterglow  phase, 
primarily  as  a  result  of  the  large  differences  in  diffusion  times 
between  the  two  chambers.  Without  an  external  source  of 
ionization,  the  much  shorter  diffusion  time  in  the  discharge 
chamber  causes  the  plasma  density  to  drop  rapidly  below  that 
in  the  ballast  volume.  This  is  now  opposite  to  the  situation  that 
exists  during  the  active  phase.  During  the  afterglow,  plasma 
from  the  ballast  volume  is  a  large  perturbation  on  the  plasma 
characteristics  in  the  discharge  chamber  and  acts  as  a  long-lived 
source  of  charged  particles  for  the  smaller  volume.  Since  this 
subject  has  not  been  well  covered  in  the  literature,  it  would  be 
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